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The basic methods of constructing the sets of mutually unbiased bases in the 
Hilbert space of an arbitrary finite dimension are reviewed and an emerg- 
ing link between them is outlined. It is shown that these methods employ 
a wide range of important mathematical concepts like, e.g., Fourier trans- 
forms, Galois fields and rings, finite and related projective geometries, and 
entanglement, to mention a few. Some applications of the theory to quantum 
information tasks are also mentioned. 
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1 INTRODUCTION 

C 

£3 . 

qh and more branches of pure mathematics, such as number theory, abstract 

algebra and projective geometry. This paper focuses on one of the most 
prominent issues in this respect, namely the construction of sets of mu- 
tually unbiased bases (MUBs) in a Hilbert space of finite dimension. For 
an updated list of open problems related to the development of quantum 
technologies the reader is directed to the Quiprocon website P . 

To begin with, we recall that two different orthonormal bases A and 
B of a d-dimensional Hilbert space TL d with metrics (. . . | . . .) are called 
mutually unbiased if and only if |(a|6)| = \j\fd for all a^A and all b^B. 
An aggregate of MUBs is a set of orthonormal bases which are pairwise 



Problems pertinent to quantum information theory are touching more 
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mutually unbiased. The MUBs have been first studied by Schwinger in 1960 
[2J. Two decades later, important results by Alltop [3] passed unnoticed 
and even a well-published paper by Ivanovic 0] still did not trigger their 
systematic research although he proved that the maximum number of such 
bases is d + 1 when d is a prime. The latter began with the important 
paper of Wootters and Fields in which they showed how to construct the 
maximum number of such bases, i.e. d + 1, for d a power of a prime. Yet, 
a still unanswered question is if there are non-prime-power values of d for 
which this bound is attained. It has been surmised [0] [I] that the maximum 
number of such bases, N(d), is equal to 1 + min(p^ i ), the latter quantity 
being the lowest factor in the prime number decomposition of d, d = Y\i p\ % 
(for a violation of this bound, see the recent work of Wocjan and Beth [H] 
and/or our comment at the end of Section 4 herein). For example, it is 
still not known [5] whether there are more than three MUBs for d = 6, the 
lowest non-prime-power dimension, although the latest findings of Wootters 
jlUj (and an earlier result of G. Tarry quoted in the last reference) seem 
to speak in favor of this conjecture. Klappenecker and Rotteler (7j showed 
that at least 3 MUBs exist in any dimension and some conditions for the 
existence of more than 3 MUBs for any dimension are also known 

MUBs have already been recognized to play an important role in quan- 
tum information theory. Their main domain of applications is the field of 
secure quantum key exchange (quantum cryptography). This is because 
any attempt by an eavesdropper to distinguish between two non-orthogonal 
quantum states shared by two remote parties will occur at the price of 
introducing a disturbance into the signal, thus revealing the attack and al- 
lowing to reject the corrupted quantum data. Until recently, most quantum 
cryptography protocols have solely relied, like the original BB84 one, upon 
1-qubit technologies, i.e., on the lowest non-trivial dimension {d = 2), usu- 
ally the polarization states of a single photon, or other schemes such as the 
sidebands of phase-modulated light [T2|. But security against eavesdropping 
has lately been found to substantially increase by using all the three bases 
of qubits, employing higher dimensional states, e.g. qudits |13|.|14j. or even 
entanglement-based protocols [T5j . Another, closely related, application of 
MUBs is the so-called quantum state tomography, which is thought to be 
the most efficient way to decipher an unknown quantum state pQ. 

Quantum state recovery and secure quantum key distribution can also be 
furnished in terms of so-called positive operator valued measures (POVMs) 
which are symmetric informationally complete (SIC-POVMs) |16j . These 
are defined as sets of d 2 normalized vectors a and b such that = 
1/Vd + 1, where a ^ b, and they are connected with MUBs. Unlike the lat- 
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ter ones, however, the SIC-POVMs are likely to exist in all finite dimensions 
and they have already been constructed for rf = 6 |]. The intricate link 
between MUBs and SIC-POVMs has recently been examined by Wootters 
jlUj and acquired an intriguing geometrical footing in the light of the "SPR 
conjecture" stating that the question of the existence of a set of d + 1 
MUBs in a d-dimensional Hilbert space if d differs from a power of a prime 
number is equivalent to the problem of whether there exist projective planes 
whose order d is not a power of a prime number. Also, Bengtsson and Er- 
icsson ^H] provided a connection with sets of d 2 facets of convex polytopes 
of power of prime dimensions: the centers of the facets can form a regular 
simplex if and only if there is an affine plane of order d that exists only if d 
is a power of a prime. 

We also mention the interesting fact recently noticed in the Lie algebra 
approach to MUBs ^Hl that a complete collection of MUBs in C d gives 
rise to a so-called orthogonal decomposition (OD) of sld{C) for which a 
longstanding conjecture says that ODs of s/^(C) can only exist if d is a prime 
power. On the other hand, for recent works on the relationship between 
MUBs and the SU(2) theory of quantum angular momentum the reader is 
directed to [2Uj |2J (in j^J a compact formula for MUBs is given in the case 
where d is a prime number). 

This survey is organized as follows. In Sections |2 and |3] the construction 
of a maximal set of MUBs in dimension d = p m , p being a prime, as a 
quantum Fourier transform acting on a Galois field {p odd) and a Galois ring 
GR(4 m ) (p = 2) is discussed. This puts in perspective the earlier formulas 
by [5| and [7], respectively. The case of non-prime-power dimensions is 
briefly examined in Section 01 Next, in Section we focus on our recent 
conjecture on the equivalence of two problems: the surmised nonexistence of 
projective planes whose order is not a power of a prime and the suspected non 
existence of a complete set of MUBs in Hilbert spaces of non-prime-power 
dimensions. The geometry of qubits is discussed and the concept of a lifted 
Fano plane is introduced. Finally, an intricate relationship between MUBs 
and maximal entanglement is emphasized in Section H3 which promises to 
shed fresh light on newly emerging concepts such as the distillation of mixed 
states and bound entanglement The entanglement properties of MUBs 
for systems of three and four qubits have been recently discussed in detail 
by Romero and collaborators We endeavored to make the paper as 

self-contained as possible from our standpoint. Yet, the interested reader 
may find it helpful to consult some introductory texts on quantum theory in 
a finite Hilbert space and its relation to Fourier transforms and phase space 
methods, e.g., the reviews by A. Vourdas |24| . 
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2 MUB'S, QUANTUM FOURIER TRANSFORMS 
AND GALOIS FIELDS 



In order to see the close connection between MUBs and Fourier trans- 
forms, we consider an orthogonal computational basis 
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with indices n in the ring 2>n of integers modulo d. There is a dual basis 
which is defined by the quantum Fourier transform 
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(2) 



where k £ Zd,Wd = exp(-^) and i = — 1. 

In the context of quantum optics this Fourier transform relates Fock states 
| A;) of light to the so-called phase states \9k}- The properties of the quantum 
phase operator underlying this construction have extensively been studied 
and found to be linked to prime number theory |25j . 



2.1 d = 2: The quantum gates approach 

For d = 2, i.e. the case of qubits, one has uj = —1 and so 

|0 O ) = -L(|O) + |1)); |^) = -L(| )-|1)). 



(3) 



These two vectors can also be obtained by applying the Hadamard matrix 
H = -j= ^ to the basis (|0),|1)). Note that the two orthogonal 

1 — 1 

bases Bo = (|0), |1)) and B\ = (|#o), are mutually unbiased. The third 
basis B2 = (\ipo), which is mutually unbiased to both Bq and B\ is 

obtained from H by the pre-action of a it/ 2 rotation S 







so that 



HS =T2 



. The three matrices (J, H, HS) thus generate the three 
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MUBs. These matrices are also important for two qubits gates in quantum 
computation |T3] . 



2.2 d = 2: The Pauli matrices approach 

The above-outlined strategy for finding MUBs for qubits contrasts with 
that used by the majority of authors. The eigenvectors of the Pauli spin 
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where = ia x a z , are precisely the sought bases Bq, B\ and Bi- 
A natural generalization of Pauli operators a x and cr^ for an arbitrary di- 
mension d is the Pauli group of shift and clock operators: 

X d \n) = \n+l), (5) 
z d\n) = u)2\n). 

For a prime dimension d = p, it can be shown that the eigenvectors of 
the unitary operators (Z p , X p , X p Z p , • • • ,X p Z p ~) generate the set of d + 
1 MUBs 26 . A natural question here emerges whether this method can 
straightforwardly be generalized to any dimension. 



2.3 MUBs on Galois fields of odd characteristic 

Let us attempt to rewrite Eq. © in such a way that the exponent of 
ojd now acts on the elements of a Galois field G = GF(p), the finite field of 
integers modulo an odd prime p (or finite field of odd prime characteristic 
p) and cardinality d = p, which in general are defined by an irreducible 
polynomial (see Sect. 3). Denoting "ffi" and "0" the two usual operations 
in the field and replacing oj^ by the root of unity u> p , we get 

|^> = -^E<»- ( 6 ) 

* n=0 

Next, we employ the Euclidean division theorem for fields |27| . which says 
that given any two polynomials k and n in G there exists a uniquely deter- 
mined pair a and b in G such that k = a n © 6, deg b < deg a. This allows 
for the exponent in Eq. 0, E, to be written as E = (a n © b) n. In 
the case of prime dimension d = p, E is an integer and the sum in Eq. Q 
is well defined. 

For G = GF(p m ), the field is defined by an irreducible polynomial, the 
cardinality is d = p m and E is a polynomial too. In this case, instead of 
Eq. © one should use the following expression 

d-l 

\ei) = _L V< [(a0ne6)0nl |n>, (7) 
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where "ir" stands for the trace of GF(p m ) down to GF(p) 



tr{E) = E@EV@ ■■■ ®E p< 



m— 1 



E G GF(p m ). 



(8) 



If p is odd, Eq. (J7J) defines the set of d bases, with the index a for the 
basis and the index b for the vector in the bases mutually unbiased to each 
other and to the computational basis Bq as well. In a slightly different 
form, this equation was first derived by Wootters and Fields 0- Its nice 
short elucidation, based on Weil sums, is due to Klappenecker and Rotteler 
[7]. Another, a more tricky derivation still in the spirit of Fourier transforms 
and claimed to hold also for the case of characteristic 2, was found by Durt 
|28| . An interesting approach based on the Weyl operators in the L 2 -space 
over Galois fields is also worth mentioning |29j . 

As already pointed out by Wootters and Fields [Sj, the reason why (J7J) 
defines the complete set of MUBs relies on the field theoretical formula 



Yln=o ut p^ a&n ^ Gn ^\ n }\ = P 1 ^ 2 ' w ith o/O andp being an odd prime. This 



method, however, fails for characteristic two where | ^n=o ^2 ' ™ l n ) I = 
for any a, b. As shown in Sect. 3 below, for this characteristic one has to 
use Galois rings instead of Galois fields to get a complete set of MUBs. 

A closer inspection of (J7J) reveals an intricate relation between MUBs 
and quantum phase operators. It is known [2Hj that the Fourier basis \9k) 
can be derived in terms of the eigenvectors of a quantum phase operator 
with eigenvalues 9k and given by 0^ = Ylt~o @k\9k)(9k\- Similarly, using 
well known properties of the field trace, one can show that each base of 
index a can be associated with a quantum phase operator 



with eigenvectors \9£) and eigenvalues 9£; the latter may thus be called an 
"MUB operator" . 



3 MUBs FOR EVEN CHARACTERISTIC FROM 
GALOIS RINGS 



Our next goal is to find a Fourier transform formulation of MUBs in 
characteristic 2. This is a very important case since the '2 m — dits' are the 
basic information units in quantum computation. One may be tempted to 
connect the Galois field algebra with the generalized Pauli operators Q by 



d-l 




(9) 



6=0 
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constructing discrete vector spaces over the Galois field |30| . For the one 
qubit case we already know that the eigenvectors of Pauli matrices a z , a x 
and a x a z define the three MUBs. Passing to the quartit (i.e., 4-dit) case, 
one finds that the operators of the following tensorial products a z a x , 
o z ® cr x a z and o x a x ® a z are associated to translations, i.e., to a single one- 
dimensional subspace in the corresponding vector space, and they define a 
unique basis represented by their simultaneous eigenvectors. Since there are 
4 + 1 one-dimensional subspaces in this discrete vector space, there are also 
4 + 1 MUBs. Other geometrically inspired derivations based on the tensorial 
decomposition of operators in the Pauli group can be found in the literature 

[HUSHES]. 

On the other hand, we are interested in a result equivalent to the for- 
mulas given in Eqs. © and for the case of characteristic two. Instead 
of the Euclidean division in the field GF(2 m ), it is necessary to consider a 
decomposition in the Galois ring GR{A m ) (defined below) so that the rele- 
vant root of unity in the Fourier formula now reads W4 = exp(2z7r/4) = i. 
For qubits GR(A) = Z4, and in general any number k in Z± can be written 
as k = a © 2 b, where © and now act in Z4. 

One now needs to introduce some abstract algebra. First one recalls that 
the Galois field GF{p m ) is the field of polynomials defined as the quotient 
Z p (x) I (q(x)) of the ring of polynomials Z p (x) by a primitive polynomial of 
order tn over Zp — GF{p). By definition, this primitive element, (^(x), has 
the property to be irreducible over the basic field GF(p), i.e., it cannot be 
factored into products of less-degree polynomials; it is also primitive over 
GF(p) of order p — 1 in the sense that it has a root a which generates any 
non zero element of GF(p) by a power sequence (a , a 2 , • • • , q p_1 = 1) and 
in addition all of its roots are in the extension field GF(p m ). There is at 
least one primitive polynomial for any extension field GF(p m ). For p = 2 
and m = 2, 3 and 4 they are, for example, of the form q(x) = x 2 + x + 1, 
x 3 + x + 1 and x 4 + x + 1, respectively. 

A Galois ring GR{A m ) of order m is a ring of polynomials which is 
an extension of Z^ of degree m containing an r-th root of unity, where 
r = 2 m — 1 t 32] [HI]- Let h,2{x) G Z%{x) be a primitive irreducible polynomial 
of degree m. There is a unique monic polynomial h(x) E Z±{x) of degree 
m such that h(x) = /i2(a:)(mod 2) and h(x)(mod 4) divides x r — 1. The 
polynomial h(x) is the basic primitive polynomial and defines the Galois 
ring GR(4: m ) = Z^(x) / (h(x)) of cardinality 4 m . This ring can be found as 
follows. Let h,2(x) = e(x) — d(x), where e(x) contains only even powers and 
d(x) only odd powers; then h(x 2 ) = ±(e 2 (x) — d 2 (x)). For m = 2, 3 and 4 
one gets h(x) = x 2 + x + 1, x 3 + 2x 2 + x — 1 and x 4 + 2x 2 — x + 1, respectively. 
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Any non zero element of GF(p m ) can be expressed in terms of a single 
primitive element. This is no longer true in Gi?(4 m ), which contains zero 
divisors. But in the latter case there exists a nonzero element £ of order 
2 m — 1 which is a root of the basic primitive polynomial h(x). Any element 
j3 G GR(4 m ) can be uniquely determined in the form (3 = a © 2 6, where 
a and b belong to the so-called Teichmiiller set T m = (0, ■•■ , £ 2m_2 ). 
Moreover, one finds that a = /3 2 " 1 . We can also define the trace to the basis 
ring i?4 by the map 

2 m -l 

trGS) = £ (10) 

fc=0 

where the summation runs over the elements of the Teichmiiller set and the 
Frobenius automorphism a reads 

ct(o©2©6) = a 2 ©2©6 2 , (11) 

with a 2 = a a. Using the 2-adic decomposition of k in the exponent of © 
and the above-given trace map, we finally get 

2 m_ 1 

\m = jM(offi206)©n]| n \. / 12 \ 



, /om ' 

the last expression gives a set of d = 2 m bases with index a for the basis and 
index b for the vectors in the basis, mutually unbiased to each other and to 
the computational base B |7j. 

Let us apply this formula to the case of quartits. In GR(4 2 ) = Z±[x\/ '(x 2 + 
x + 1) the Teichmiiller set reads T 2 = (0, 1, x, 3+3x); the 16 elements a©2©6 
with a and b in I2 are shown in the following matrix 



2 2x 2 + 2x 

1 3 l + 2x 3 + 2x 
x 2 + x 3x 2 + 3a; 

3 + 3x 1 + 3x 3 + x 1 + x 



(13) 



Extracting the Teichmiiller decomposition (a © 2 b) n = a' © 2 b' and 
calculating the exponent tr(a' © 2 b') = a' © 2 b' © a' 2 © 2 6' 2 one gets 
the four MUBs 

Bt = (1/2){(1, 1, 1, 1), (1, 1, -1, -1), (1, -1, -1, 1), (1, -1, 1, -1)} 
B 2 = (1/2){(1, -1, -%, -i), (1, -1, i, i), (1, 1, i, -i), (1, 1, -i, *)} 
B 3 = (1/2){(1, -*, -i, -1), (1, -i, i, 1), (1, i, i, -1), (1, i, -i, 1)} 
fi 4 = (1/2){(1, -i, -1, -*), (1, -i, (1, 1> -0, (1, -1, 01- (14) 
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The case of 8-dits can be examined in a similar fashion, with the ring 
GR(4 3 ) = Zi[x]/(x 3 + 2x 2 + x — 1) and Teichmiiller set featuring the fol- 
lowing eight elements: T2 = {0, 1, x, . x 2 , 1 + 3x + 2x 2 , 2 + 3x + 3x 2 , 3 + 3x + 
x 2 ,l + 2x + x 2 }. 



4 MUB'S FOR NON-PRIME-POWER DIMEN- 
SIONS 

Here we shall discuss in some detail the simplest case, which is d = 6, 
the lowest non-prime-power (n-p-p) dimension. For this case, one constructs 
a set of three MUBs as follows. One takes the three MUBs in d = 2, viz. 



= (|0>, |1», b« = (|#o>, B a w = (IVo), IVi}), 



(15) 



or, in the matrix form, Bq=I2,B±=H and B^ —HS, and the first three 
MUBs in d = 3, viz. 



I?, 



(2) 



1), |2», B{ 2) = Quo), |u 2 », 1^ = (|«o), h», (16) 



(2) 



or, in a more convenient form 



B, 



(2) 



(2) 



1 1 1 

1 w 3 cJ 3 
1 w 3 u 3 



(2) 



1 

7s 



1 w 3 w 3 

1 cJ 3 1 

1 1 W"3 



(17a) 



and extracts the expressions for three MUBs in d = 6 from the rows of the 



following tensorial product matrices Co = B^ (g) B^ = Iq, C\ = B| ±; ® B 1 
and C2 = B^p (g) B^ 2 "*. This construction can easily be generalized to any 
n-p-p dimension [Jj One considers the prime number decomposition 

d = nT—iPf*, takes its smallest factor m = miniipf 1 ), and gets m + 1 MUBs 
from the tensorial product B^ = ®\ = iB^\ (k = 0, ..,m). 

At this point, it is instructive to enlighten the above-described construc- 
tion of MUBs by confining ourselves to the Galois ring in d = 6. Let us 
take the latter as the quotient GR(6 2 ) = Zq[x]/ (x 2 + 3x + 1) of polynomials 
over Zq by a polynomial irreducible over both Z2 and Z3. GR(6 2 ) has 36 
elements. The notion of Teichmiiller set can be generalized to the so-called 
Sylow decomposition [£]. Any element G GR(Q) can be uniquely deter- 
mined in the form (5 = a b, where a and b are in the Sylow subgroups 
S a and Sb- These can be defined as S a = {x G GR(6) : 2x = 0} and 
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S b = { x G GR{6) : 3x = 0}, i.e. 



S a = {0, 3,3x,3 + 3x}, 
S b = {0, 2, 4, 2x, Ax, 2 + 2x, 2 + 4x, 4 + 2x, 4 + 4x}. (17) 

Since the quotient polynomial is irreducible, one notices that S a and Sb 
themselves are finite fields, being isomorphic to GF(4) and GF(9), respec- 
tively. One can therefore express the ring in dimension 6 as GF(4)(BGF(9) = 
GR(6). Can this property be useful to construct MUBs themselves, or it 
merely represents a constraint on the maximum number of MUBs? One 
construction of MUBs for d=6 was based on the tensorial product of MUBs 
in dimension 2 and 3, respectively. But the three MUBs in dimension two 
do not follow from the four elements of GF(4), but from the four elements 
of Gi?(4 x ) = Z±. On the other hand, the four MUBs in d=3 follow from the 
three elements of GF(3) = Z3. So the decomposition of GR(6) as a prod- 
uct of two fields appears to be irrelevant to the topic of MUBs. Moreover, 
it was shown that complete sets of MUBs in n-p-p dimensions cannot be 
constructed using a majority of generalizations of known formulas for finite 
rings [Hj. This, however, should not deter us from looking at other pos- 
sible constructions. For example, using the properties of sets of mutually 
orthogonal Latin squares, it has recently been shown that in the particular 
square dimension d = 26 2 it is, in principle, possible to construct at least 
6 MUBs, while the construction based on the prime number decomposition 
determines only mini (p? j ) + l = 2 2 + l = 5of them [S] . 

5 MUB'S AND FINITE PROJECTIVE PLANES 

An intriguing similarity between mutually unbiased measurements and 
finite projective geometry has recently been noticed Let us find the 

minimum number of different measurements we need to determine uniquely 
the state of an ensemble of identical d-state particles. The density ma- 
trix of such an ensemble, being Hermitic and of unit trace, is specified by 
(2d 2 /2) — 1 = d 2 — 1 real parameters. When one performs a non-degenerate 
orthogonal measurement on each of many copies of such a system one even- 
tually obtains d — 1 real numbers (the probabilities of all but one of the d 
possible outcomes). The minimum number of different measurements needed 
to determine the state uniquely is thus (d 2 — l)/(d — 1) = d+ 1 [30] . 

It is striking that the identical expression can be found within the con- 
text of finite projective geometry. A finite projective plane is an incidence 
structure consisting of points and lines such that any two points lie on just 
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one line, any two lines pass through just one point, and there exist four 
points, no three of them on a line From these properties it readily 

follows that for any finite projective plane there exists an integer d with 
the properties that any line contains exactly d + 1 points, any point is the 
intersection of exactly d + 1 lines, and the number of points is the same as 
the number of lines, namely d 2 + d + 1. This integer d is called the order 
of the projective plane. The most striking issue here is that the order of 
known finite projective planes is a power of prime. The question of which 
other integers occur as orders of finite projective planes remains one of the 
most challenging problems of contemporary mathematics. The only "no-go" 
theorem known so far in this respect is the Bruck-Ryser theorem |37| saying 
that there is no projective plane of order d if d — 1 or d — 2 is divisible by 4 
and d is not the sum of two squares. Out of the first few non-prime-power 
numbers, this theorem rules out finite projective planes of order 6, 14, 21, 
22, 30 and 33. Moreover, using massive computer calculations, it was proved 
that there is no projective plane of order ten. It is surmised that the order 
of any projective plane is a power of a prime. 

It has been conjectured by Saniga and two of us |T2] that the question 
of the existence of a set of d + 1 MUBs in a d-dimensional Hilbert space if 
d differs from a power of a prime number is identical with the problem of 
whether there exist projective planes whose order d is not a power of a prime 
number. Furthermore, for power of a prime dimension, the complete sets 
of MUBs can be put in correspondence with d + 1-arcs, which are 'curves' 
known as ovals in (Desarguesian) projective plane of order d |38j . For d = 2 n 
and n > 3 there are two types of ovals, viz. conies and non-conics, implying 
the existence of two types of MUBs for such dimensions. In addition, in the 
same case of a power of a prime dimension d = p r , the p r vectors of a basis 
correspond to the total number of points in a so-called neighbour class along 
a (proper) conic of a projective Hjelmslev plane defined over a Galois ring 
of characteristic p 2 and rank r, whereas the d+1 MUBs correspond exactly 
to the total number of pairwise disjoint neighbour classes on the conic |39| . 

5.1 GF(8) and the Fano plane 

The smallest projective plane, also called the Fano plane, is obviously 
the d = 2 one; it contains 7 points and 7 lines, any line contains 3 points 
and each point is on 3 lines. It comprises a 3-dimensional vector space over 
the field GF(2), each point being a triple (<7i, <?2>53)> excluding the (0,0,0) 
one, where gi G GF(2) = {0, 1} |36j . The points of this plane can also be 
represented in terms of the non-zero elements of the Galois field G = GF(2 3 ). 
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To see this, we recall that this field is isomorphic to Z<i (x)/(p(x)) with 
the polynomial p(x) = x 3 + x + 1 irreducible in GF{2). It is well-known that 
there are three useful representations of the elements of GF(8) as shown in 
Table 1 [SHI HOI HU 

Table 1: Representations of the elements of the Galois field GF(8) 



as powers of a 


as polynomials 


as 3-tuples in Z\ 








(0,0,0) 


1 


1 


(0,0,1) 


a 


X 


(0,1,0) 


a 2 


X 


(1,0,0) 


a A 


1 + x 


(0,1,1) 


a 4 


1) 

X + X 


(1,1,0) 


a b 


1+x+x 2 


(1,1,1) 


a 6 


1 + x 2 


(1,0,1) 



The first representation emphasizes the fact that G* = G — {0} is a 
multiplicative cyclic group of order 7, for a 7 = 1 . The second representation 
is obtained from the first by calculating modulo the primitive polynomial 
p{x). Finally, the 3-tuple representation is obtained from the coefficients of 
the three powers x° = 1, x 1 = x and x 2 . Taking these 3-tuples as the points 
of a 3-dimensional vector space, we recover the Fano plane — see Fig. 1. 

5.2 The lifted Fano plane in GR(A 3 ) 

We already know from Sect. 01 that the relevant object for 2 m -dits is not 
the Galois field GF(2 m ), but rather the Galois ring GR(4 m ). It is therefore 
important to have a look at the geometry in the space A = GR(4 3 ). For a 
ring, the concept of a vector space must be replaced by that of a module. The 
largest cycle in A is the set T£ = T% — {0} (see Sect. 3), and each element of 
Xj* can be represented in the same way as in the case of a Galois field. This 
is summarized in Table 2. Any polynomial h(x) in (column 2) is uniquely 
projected as a polynomial h,2{x) = h(x) (mod 2) in GF(8), which results in 
the 3-tuple representation in Z\ (column 4). Vice versa, any polynomial in 
GF(8) has a unique lift in . Since the geometrical structure we are looking 
at is combinatorial and doesn't depend on particular coordinates, it follows 
that the lifted Fano plane in T£ is still the Fano plane up to isomorphism. 
So the Fano geometry is inherent in the geometry of qubits, but we needed 
a special coordinatization in order to be able to see that. The very recent 
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Table 2: Representations of the elements of the cyclic group in the Galois 
ring GR(4 3 ) 



as powers of £ 


as polynomials 


as 3-tuples in Zf 


as 3-tuples in Z\ 








(0,0,0) 


(0,0,0) 


1 


1 


(0,0,1) 


(0,0,1) 




X 


(0,1,0) 


(0,1,0) 


e 


o 

X 


(1,0,0) 


(1,0,0) 


e 


1 + 3x + 2x 2 


(2,3,1) 


(0,1,1) 


£ 4 


2 + 3x + 3x 2 


(3,3,2) 


(1,1,0) 


e 


3 + 3x+x 2 


(1,3,3) 


(1,1,1) 


e 


1 + 2x + x 2 


(1,2,1) 


(1,0,1) 



work of Planat, Saniga, and Kibler |42j relates the understanding of MUBs, 
the Fano plane, entanglement, and quantum paradoxes to the construction 
of projective lines over special families of finite rings. 



(1,1,0) 




(1,0,0) (0,0,1) (1,0,1) 



Figure 1: The Fano plane. 
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6 ORTHOGONAL SETS OF MAXIMALLY 
ENTANGLED STATES 



The above-discussed methods of constructing MUBs can straightfor- 
wardly be used for recognizing orthogonal partial bases of maximally en- 
tangled states, of which some can be mutually unbiased in their correspond- 
ing subspaces implying that the MUB normalization 1/yd is replaced by a 
higher normalization l/y/cL, with d s < d. In this section, by partial bases we 
mean orthogonal sets of Hilbert vectors that are not in enough number to 
form complete bases. Following the methodology outlined in Sections El and 
01 let us consider a set of generalized Bell states defined as a two particle 
quantum Fourier transform ^1] |43| 

d-l 

\B h ,k) = -^Y, u d n \n,n + h), (18) 

n=0 

where \n, n+h) denotes the two-particle state \n)\n+h) and the operation n+ 
h is performed modulo d. These states are both orthonormal, (Bh,k\Bh' ,k') = 
Shh'Skk', and maximally entangled, trace2\Bh,k)(Bh,k\ = \ld, where trace2 
means the partial trace over the second qudit |13j . If one restricts to the 
case of 2-qubits, one recovers the well-known representation of Bell states 
(|Bb,o>, |B 0)1 » = ^(|00) + |11), |00> - |11)), (|Bi, >, |Bi,i»= 
^(|01) + |10),|01) - 1 10>), where a more compact notation |00) = |0,0), 
|01) = |0, 1),. . . , is employed. Let us first focus on the 2-qubit case. Paral- 
leling what we did in Section |31 kn in l|18|) is first identified as the multipli- 
cation kQn of polynomials in GR(4) and then k is Teichmiiller decomposed, 
i.e., k = a © 2 b. This leads to 4 sets (h, a = 0,1) of two vectors (b = 0,1), 
namely 

\Bl b ) = ^jz^ ameb)Qn \n,n@h). (19) 

* n=0 

Casting the last equation into its matrix form (save for the proportionality 
factor), 



(|00» + |11), |00> - |11)); (|01) + 1 10), |01> - |10» 
(|00) + |00) - i|ll)); (|01) + »|10), |01>) - i|10)) 



(20) 



one finds that two partial bases in one column are mutually unbiased in 
their subspace, while the vectors in two partial bases on the same line are 
orthogonal to each other. 
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Eq. H19|) can easily be extended to maximally entangled two-particle sets 
of 2 m -dits by applying, as in Eq. ((T2*|) . the Frobenius map ifTUj) to the basis 
field i?4 

\B a h , h ) = ^= 2 ^^ 2Gb ^\n,n®h). (21) 

For 2-particle sets of quartits, using Eqs. ((Tl)) and (|2*Tj) . one thus gets 4 sets 
(|£>^ b ), h = 0, 3) of 4 mutually unbiased partial bases (a = 0, 3), 

{(|00) + |11) + |22) + |33), |00) + |11) - |22) - |33), 
|00) - |11) - |22) + |33), |00) - |11) + |22) - |33»; 
(|00) - |H) - i\22) - i|33>, |00> - |11) + i\22) + »|33), 
|00) + |H) + i\22) - i\33), |00) + |11) - i\22) + «|33»; 

•••} 

{(|01> + |12) + |23) + |30), |01) + |12) - |23) - |30), 
|01) - |12) - |23) + |30), |01) - |12) + 1 23) - |30>); 
(|01) — 1 12> — »|23) - i|30), |01) - 1 12) + i\23) + i|30), 
|01) + 1 12) + i\23) - i\30), |01) + \12) - i\23) + «]30>); 

•••} 

{(|02) + 1 13) + |20) + |31), |02) + 1 13) - |20) - |31), 
]02) — |13> — |20) + |31), |02) - 1 13) + |20) - |31>); • • • 

•••} 

{(|03) + 1 10) + |21) + |32), |03) + 1 10) - |21) - |32), 
|03) — |10> — |21) + |32), |03) - 1 10) + |21) - |32)); • • • 

•••}, (22) 

where we have skipped the partial normalization factor (1/2). Within each 
set, the four partial bases are mutually unbiased, as in (|14j) . while the vectors 
of the partial bases from different sets are orthogonal. 

Turning now to odd characteristic, i.e. to d = p m with p an odd prime, 
we can similarly extend Wootters formula (jJJ) to the generalized Bell states 

W = ^E"T &nm&n] \^®h), (23) 
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where the trace is defined by Eq. (jSJ) . A list of the generalized Bell states of 
qutrits for a = can be found in |44| . a work that relies on a coherent state 
formulation of entanglement. In general, for d a power of a prime, starting 
from (|18[) or (]23|) one obtains d 2 sets of d maximally entangled states. Each 
set of the d bases (with h fixed) has the property of mutual unbiasedness. 

Eq. ((TH|) can be used, without any substantial restriction, to find d sets 
{h = 0, ..,d — 1) of maximally entangled states in any composite dimen- 
sion d = YYi=iPi ■ Or one can also follow the strategy of Section 0] to get 
777, = mini^') sets of MUBs of maximally entangled states. In d = 6, for 
example, one expects that two such d sets can be constructed. Using the 
tensorial products in Sect. 0J one indeed finds the two 2x6 sets (with the 
1 / y/G factor omitted) 

{(|00) + |11) + |22) + |33) + |44) + |55>, |00) + w 3 |ll) 4- w 3 |22) + |33> + u; 3 |44) + w 3 |S5), 
|00> +" 3 |H> + w 3 |22) + |33) + w 3 |44) + u 3 |55), |00> + |11> + |22> - |33> - |44) - |55>, 
|00> + w 3 |ll> + w 3 |22) - |33> - u> 3 [44) - l7 3 |55>, |00> +w 3 |ll> + w 3 |22) - |33) - w 3 |44) - w 3 |55»; 
(|00) + w 3 |ll) + u 3 |22) + i|33) + i^ 3 |44) + iuj 3 \55), (00) + oJ 3 |ll> + |22> + i|33) + ioJ 3 |44> + i|55>, 
|00> + |11> + oJ 3 |22) + s|33) + i|44) + iu, 3 |55), 100) +w 3 |ll> + w 3 |22) - *|33) - ia> 3 |44) - iw 3 |S5), 
|00} +(J- 3 |11) + |22> - i|33) - iw 3 |44) - i|55), |00> + |11> + w 3 |22) - i|33> - z|44> - »cJ* 3 |55»i . . .} 



{(|01> + 12) + 23) + [34) + |4S) + |50>, |01) + w 3 |12) + w 3 |23> + |34> + ^ 3 |45) 4-cr 3 |50>, • • ■ }. (24) 

Multipartite entanglement is a key ingredient of many quantum proto- 
cols, still needing much work to be properly understood. Sets of orthogonal 
product states that are not extendible, meaning that no further product 
states can be found orthogonal to all the existing ones, have recently at- 
tracted a lot of attention. These non-extendible product bases [45] . and their 
complement [46] . certainly deserve reconsideration in terms of the above- 
outlined theory, which is based on abstract algebra and finite geometry. 

The Fourier transform approach implies that mutual unbiasedness and 
maximal entanglement are complementary aspects in orthogonal quantum 
measurements. In such measurements, the quantum states are encoded in a 
three-dimensional lattice of indices h (entanglement), a (unbiasedness) and 
b (dimensionality of Hilbert space). If d is a power of a prime, the lattice is 
a cube since in this case h, a and b reach their limiting value d. If one for- 
gets about entanglement (h = 0), the finite geometry which seems to be of 
most relevance is that of a finite projective plane. On the other hand, when 
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unbiasedness is not taken into account, as well as for multipartite informa- 
tion tasks when d is not (a power of) a prime, other concepts have been 
introduced, such as Bell inequalities [37], coherent states [13], entanglement 
swapping 02] > generalized Hopf fibrations jlS] , topological entanglement |50| 
and bound entanglement |45| . to mention a few. 
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